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This investigation evaluates the numerical prediction of flow distortion and pressure recovery for a bound- 
ary layer ingesting offset inlet with active flow control devices. The numerical simulations are computed using 
a Reynolds averaged Navier-Stokes code developed at NASA. The numerical results are validated by compar- 
ison to experimental wind tunnel tests conducted at NASA Langley Research Center at both low and high 
Mach numbers. Baseline comparisons showed good agreement between numerical and experimental results. 
Numerical simulations for the inlet with passive and active flow control also showed good agreement at low 
Mach numbers where experimental data has already been acquired. Numerical simulations of the inlet at high 
Mach numbers with flow control jets showed an improvement of the flow distortion. Studies on the location 
of the jet actuators, for the high Mach number case, were conducted to provide guidance for the design of a 
future experimental wind tunnel test. 


I. Nomenclature 

AIP Aerodynamic Interface Plane 

Aoo Inlet capture area, (pooKoAoo = p aipVaipAaip) 

Aaip Area of inlet at AIP 

BL Boundary Layer 

BLI Boundary Layer Ingesting (or Ingestion) 

BWB Blended Wing Body 

CFD Computational Fluid Dynamics 

Cm j e t Total jet momentum coefficient, Cm jet = (total jet momentum) /(poo Aaip ) 

D AIP diameter 

DC60 Distortion using a 60° sector 

DPCP avg Average SAE circumferential distortion descriptor 

L Inlet length 

M Mach number 

P Pressure 

Pp Total Pressure 

Ppavg Average Total Pressure at the AIP 

(. Pravg)crit Average Total Pressure of the ’worst’ 60° sector 

h Vane height 

MFR Mass Flow Ratio, (Total jet actuator mass flow rate )/( Inlet mass flow rate) 

PAVi Average total pressure of ring i 

PAV LOWi Average low total pressure of ring i 
q avg Mean dynamic pressure at the AIP 
Rep> Reynolds number based on AIP diameter D 
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scfm 

Standard Cubic Feet per Minute 

Uoo 

Free- stream Velocity, m/s 

u,v,w 

Velocity in x,y,z directions respectively, m/* 

VG 

Vortex Generator 

x,y,z 

Cartesian axes 

5 

Boundary layer height 

Poo 

Free- stream Density 


A. Subscripts 

avg Average 

oo Free- stream 


II. Introduction 

I N an effort to reduce the environmental impact of commercial aircraft using revolutionary propulsion technologies, 
NASA initiated the Ultra Efficient Engine Technology (UEET) program. One of the elements of the UEET program 
is the application of flush mounted, boundary layer ingesting (BLI), offset inlets on the aft portion of an aircraft. 
System studies for the Blended Wing Body (BWB) transport have shown significant reductions in fuel burn by using 
this type of inlet. For the BWB vehicle, a BLI inlet placed on the upper rear surface of the wing would have a 
boundary layer to inlet height ratio of 30%. The ingestion of such a large boundary layer coupled with the S-shaped 
offset of the inlet diffuser, results in a large flow distortion at the engine fan face. ’ ’ Experiments have shown that 
inlet distortion can be improved for the ingestion of a 30% thick boundary layer to acceptable levels using flow control 
devices located inside the inlet. 

The current investigation studies the effects of vortex generating (VG) vanes and fluidic jets on inlet distortion and 
pressure recovery. The objectives of this investigation are twofold. The first objective is to validate the flow solver 
for the BLI offset inlet flow with and without flow control by fully modeling the VG vanes and jets. Once validated, 
simulations of the BLI inlet, with VG jets, will be performed for a future high Mach number experiment. The results 
from these simulations will provide guidance on the placement of jets in the BLI inlet. The second objective of this 
study is to provide additional insight into the flow physics of VG vanes and jets inside a BLI offset inlet. The insights 
gained from the simulations of the fully modeled control devices, will be used for the future validation of a source 
term actuator model. The source term model will significantly reduce the cost of simulating flow control devices in 
an inlet. Thus, making optimization studies for the placement and sizing of actuators in an inlet viable. 

III. Numerical Modeling 

The steady-state flow field for the BLI offset inlet with VG 
vanes and jets was computed using the flow solver code, OVER- 
FLOW, 7, developed at NASA. This code solves the compressible 
Reynolds averaged Navier- Stokes (RANS) equations using the di- 
agonal scheme of Pulliam and Chaus see. The RANS equations 
are solved on structured grids using the overset grid framework of 
Steger et al. This overset grid framework allows for the use of 
structured grids for problems which have complex geometries. To 
improve the convergence of the steady-state solution, the OVER- 
FLOW code also includes a low-Mach number preconditioning op- 
tion and a multigrid acceleration routine, which were both used for 
the numerical simulations. All of the simulations in this study used 
Menter’s two-equations (k-co) Shear-Stress Transport (SST) turbu- 
lence model. The SST turbulence model was found to be the 
best turbulence model option in OVERFLOW for the simulation 
of streamwise vortices embedded in a boundary layer. 

The numerical simulations were performed using the parallel 
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Figure 2. Close-up view of the overset nozzle grids for a YG jet Figure 3. View of the overset nozzle grids and inlet grids, 
inside of the BLI inlet. 


version of the OVERFLOW code developed by Buning. This code uses the Message-Passing Interface (MPI) and 
can run on a tightly-coupled parallel machine or a network of workstations. The code distributes zones to individual 
processors and can split larger individual zones across multiple processors using a domain decomposition approach. 

The structured overset grid system was generated using the Chimera Grid Tools package. Figure 1 shows a close- 
up view of the overset grids near the VG vanes on the inlet surface. The vanes were modeled as rectangular fins which 
was shown to be comparable to a fully modeled trapezoidal vane. Figure 2 shows a close-up view of the nozzle grid 
system for the VG jet simulation. The steady jet is skewed 90° to the frees-tream flow and pitched at an inclined angle 
of 30° to the surface. These pitch and skew angles for the jet result in the generation of a single stream wise vortex. 
This jet is simulated by modeling the nozzle plenum below the surface of the inlet. This simplifies the inflow boundary 
condition for the jet by letting the flow develop in the nozzle plenum and exiting at the duct surface. Figure 3 shows 
the inlet grids with the VG jet grids on the bottom surface of the inlet. 


IV. Wind Tunnel Experiments 

The numerical simulations of the BLI offset inlet were compared to two experimental tests conducted at NASA 
Langley Research Center. ’ One experiment was conducted in a low speed wind tunnel, the other in a high speed 
transonic tunnel. The low speed experiment evaluated both passive and active flow control devices; the high speed 
experiment investigated the baseline inlet at flight Mach and Reynolds numbers. 

A. Low Mach and Low Reynolds Number Experiment 

The low speed experiments for the BLI offset inlet were performed at NASA Langley’s Basic Aerodynamics Research 
Tunnel (BART). This wind tunnel had a free-stream Mach number of 0.15 as was able to generate a boundary layer 
to inlet height ratio of 36%. The experimental test had a Reynolds number of Rep = 0.54 • 10 6 , where D = 6.0 inches 
and is equal to the aerodynamic interface plane (AIP) diameter. This experiment showed the improvement of flow 
distortion using VG vanes and jets located inside the inlet. The data from this experiment will be used to compare and 
validate the flow solver for low Mach number flows. Further information on the experiment is provided by Gorton et 
al. 3 4 

B. High Mach and High Reynolds Number Experiment 

The high Mach and Reynolds number experiments were conducted at NASA Langley’s 0.3 -Meter Transonic Cryogenic 
Tunnel (0.3-Meter Tunnel) for the BLI offset inlet. Experimental data was obtained for the baseline inlet case for a 

3 of 18 


American Institute of Aeronautics and Astronautics Paper 2004-2318 



free-stream Mach number range from 0.25 to 0.83. The experimental test had a Reynolds number range of Reo = 6.8 • 
10 6 to 14.3 • 10 6 , where the AIP diameter, D = 2.448 inches. This experiment was able to test the BLI inlet at 
actual flight Mach and Reynolds numbers as expected for the BWB aircraft application. This experiment generated a 
boundary layer of approximately 30% of the inlet height ratio. 


C. BLI Offset Inlet Geometry 


BLI Inlet Side View 


The BLI offset inlet was designed by Boeing under 
contract to NASA in order to provide an inlet that 
would be representative of designs considered for 
the commercial version of the BWB aircraft. The 
inlet geometry is of a generic nature and is an open 
geometry that can be used for computational val- 
idation and experimental investigations. Figure 4 
shows the side view of the S-shaped inlet for the 
BART experiment. The inlet for the BART experi- 
ment is 6.129% of full scale and the inlet model for 
the 0.3 -Meter Tunnel experiment was built to 2.5% 
scale. The inlet cross section transitions from a D- 
shape at the entrance to circular shape at the AIR 
Figure 4 also shows the location of the AIP where 
the total pressure rake was located and where the 
engine face would be when attached to the inlet. 

This figure also shows the location of the VG plane 
where the passive and active flow control devices 

were placed for the BART experiment. This location has 32 jet orifices placed circumferentially on the inlet surface 
where different combinations of jet locations were evaluated by Gorton et al. 



Figure 4. BLI offset inlet geometry for the 6.129% scale BART model, 
showing the location of flow control devices and the AIP. 


D. Distortion Descriptors 

Two different distortion descriptors were used to evaluate inlet performance. The first inlet performance parameter is 
the average SAE circumferential distortion descriptor defined in the Aerospace Recommended Practice (ARP) 1420 
standard. In this paper the average distortion intensity is given by DPCP avg and defined in (1). 


DPCPavg — ^ 


Y A rings 


rings i= i 


^ Intensityi 


( 1 ) 


where i is the ring number on the AIP rake and N r i ngs , the total number of rings. The Intensity for each ring is defined 
in (2). 

, PAVt-PAVLOWi 

Intensityi = (2) 

where PAVi is the average total pressure of ring i and PAVLOWi , the average of the low total pressure region below 
PAV t . 

The second engine face distortion descriptor used is DC6 0 defined in (3). 


DC60 = ^ Tavg ^ Tav s) crit Q) 

Qavg 

The DC6 0 distortion descriptor is computed from the mean total pressure at the AIP, Pravg , the mean dynamic pressure, 
q avg , and the mean total pressure in the ’worst’ 60° sector, ( Pravg)crit - 15 Where the ’worst’ 60° sector is the sector 
with the lowest mean total pressure. Unlike the DPCP avg distortion descriptor, DC 60 is scaled by the average dynamic 
pressure at the AIP. For experimental data, q avg is typically estimated using static pressure measurements on the surface 
of the inlet around the rake. The mean of all the static pressure measurements is then used to estimate the Mach number 
at the total pressure rake locations which is then used to compute q avg . For the high Mach number cases with BLI, 
it was discovered using data from the numerical simulations, that this approach to estimating q avg resulted in a 7 to 
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Figure 5. A comparison of the numerical and experimental results for the BLI inlet for the baseline low speed case, Moo = 0.15, with a duct 
flow rate of 1800 scfm. These contour plots show the total pressure ratio at the AIP where the CFD results are interpolated onto the 120 
probe locations used in the experiment. 


14% error. This estimation error in q avg is then carried over to the DC 60 distortion values. For the low Mach number 
experiments, the estimation error for q avg was 1 to 4% since the static pressure did not vary as much across the AIP as 
in the high Mach number case. All of the DC6 0 results for the numerical simulations and experimental data followed 
the practice of computing DC6 0, where q avg is estimated from a mean static wall pressure and the total pressures at the 
rake locations. The static and total pressures used to compute the distortion descriptors for the numerical simulations 
were also interpolated at same locations as the experimental data. This way the numerical results can be evaluated at 
the same resolutions as the experimental data. 

V. Results and Discussions 

A. Baseline: Low Mach and Low Reynolds Number Case 

The baseline case for the low Mach and Reynolds number experiment had a free- stream Mach number of 0.15 with a 
Reynolds number of Reo = 0.54 • 10 5 6 . The duct mass flow rate was held fixed at 1800 scfm for all low Mach number 
experimental data presented in this paper. Figure 5 shows a comparison between the numerical and experimental 
results for the low Mach and Reynolds number baseline case. This figure shows the contour plots of the total pressure 
ratio, Pt / (Pt) «>, at the AIP. The experimental data was taken using a 40 probe total pressure rake located at the AIP. 
This rake was then rotated by 15 and 30 degrees resulting in 120 total pressure measurements. The high resolution of 
the numerical results in Fig. 5 are interpolated to the same resolution as the experimental data. The 120 total pressure 
measurements are then used to compute a distortion value using the DC60 method. The experiment had a DC60 
value of 0.30 where the numerical simulation was predicting 0.26. This comparison shows a difference of 13% for 
the low speed baseline results. Overall the contour plots compare very well for the baseline case showing similar flow 
distortion patterns. 

B. Baseline: Low Mach and High Reynolds Number Case 

Numerical simulations were compared to the high speed experiments taken in the 0.3 -Meter Transonic Cryogenic 
Tunnel at two free-stream Mach numbers. The lower free-stream Mach cases are shown in Fig. 6 where Moo = 0.25 
with mass flow rates of 6.38 and 7.10 Ibm/s. These duct flow rates result in an acceleration of the flow inside the duct 
for an average Mach of 0.35 at the AIP for the low mass flow rate and 0.39 for the high mass flow rate case. These 
two cases also have Reynolds numbers, based on the AIP diameter, of about Reo = 6.9 • 10 6 . 

Since the BLI inlet was modeled on a flat plate, the effects of the tunnel walls were not accounted for in the 
numerical simulations. In order to match the experimental flow conditions, the flat plate length ahead of the inlet and 
the free-stream Mach number were adjusted, matching the boundary layer velocity measured near of the inlet face. 
The BL rake was located 3.784 inches from the centerline of the inlet and 0.10 inches upstream of the cowl highlight. 
The rake was approximately 1.67 inches away from the cowl outer surface and had a height of 0.58 inches. In the first 
numerical simulation the flat plate length was determined to match the experimental BL height. This simulation used 
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Figure 6. Comparison of numerical and experimental results for the BLI inlet for the baseline flow taken in the 0.3-Meter Transonic 
Cryogenic Tunnel. Case (a) has a free-stream Mach of 0.25 with a Rep = 6.8 • 10 6 and a duct mass flow rate of 6.38 lbm/s. Case (b) has a 
free-stream Mach of 0.25 with a Reo = 6.9 10 6 and a duct mass flow rate of 7.10 lbm/s. 


the free-stream Mach number given by the experiment which was measured upstream of the test section. Figure 7 
shows the BL rake data for the high and low Mach number cases at a given inlet mass flow rate. These figures show 
how the Mach number at the BL edge for the numerical simulations was slightly higher than the experiment. The 
free-stream Mach number for the simulations was then adjusted to match the velocity measured at the BL rake. The 
free-stream Mach number for the Moo = 0.25 case was adjusted to Moo = 0.234 in the numerical simulation, producing 
a better match to the BL velocity as shown in Fig. 7. In the Moo = 0.833 case, the free-stream Mach number was 
reduced to Moo = 0.784 which resulted in a better match to the BL velocity. 

The BL comparison in Fig. 7 shows how the BL profile is slightly different in the experiment for the high Mach 
number case as compared to the numerical simulations. The BL in the experiment has less energy near the wall than 
the numerical simulation. The BL for the experiment was generated from the tunnel wall and not from a splitter plate, 
which may account for the difference in the BL profiles. 

The flow distortion for the high Reynolds number experiments was calculated using the DPCP avg descriptor. The 
flow distortion for the lower mass flow rate case had a DPCP avg = 0.005 as compared to 0.006 predicted by the 
numerical simulation. The numerical simulation had a slightly lower pressure recovery of Pravg/ ( Pt )oo = 0.991 which 
compares well to the experimental value of 0.994. The distortion for the higher mass flow rate was identical for the 
simulation and the experiment having at value of DPCP avg = 0.006. As in the lower mass flow rate case, the numerical 
results showed a slightly lower pressure recovery of 0.989 as compared to the experiment which measured 0.994. 
Figure 6 shows a comparison of the total pressure ratio contours which compare very well between the experiment 
and numerical results. 

C. Baseline: High Mach and High Reynolds Number Case 

The high free-stream Mach number cases are shown in Fig. 8 where Moo = 0.833 and the duct mass flow rates were 5.0 
and 6.0 lbm/s. The low mass flow case had a Rep = 14.3 • 10 6 and the high mass flow rate case a Ren = 13.8 • 10 6 . The 
comparison of distortion patterns for the high free-stream Mach number cases are good overall as seen in Fig. 8. For 
the low mass flow rate case in Fig. 8 a, a DPCP avg of 0.040 was predicted by the numerical simulations as compared 
to 0.034 measured in the experiment. The numerical simulations predict a distortion about 18% larger than was seen 
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Figure 7. A comparison of the boundary layer profile on the side of the inlet for the experiment and the numerical simulation. The rake is 
located 0.10 inches upstream of the highlight on the inlet cowl, 3.784 inches from the inlet centerline and approximately 1.67 inches from 
the outer surface of the inlet cowl. These figures show how the free-stream Mach number for the simulation needed to be adjusted to match 
the velocity profile at the BL rake. 


(a) 


(b) 


Baseline 


M =0.83 


DPCP ave =0.040 


D PC P AVE =0.034 


Pt/(Pt) 



Baseline 


M =0.83 


DPCP ave =0.057 


DPCP. 


=0.048 


PtAPtL 



Figure 8. Comparison of numerical and experimental results for the BLI inlet for the baseline flow taken in the 0.3-Meter Transonic 
Cryogenic Tunnel. Case (a) has a free-stream Mach of 0.833 with a Reo = 14.3 • 10 6 and a duct mass flow rate of 5.00 lbm/s. Case (b) has a 
free-stream Mach of 0.833 with a Reo = 13.8 10 6 and a duct mass flow rate of 6.01 lbm/s. 
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Figure 9. A comparison of the numerical and experimental results for the BLI inlet using VG vanes for the low speed case. These contour 
plots show the total pressure ratio at the AIR The CFD results were interpolated onto the 120 probe locations used in the experiment. 


in the experiment. This higher distortion value can be explained by looking at the total pressure contour plots in 
Fig. 8 a. Comparing the two plots it can be seen that the numerical simulation has a higher peak total pressure and a 
lower minimum. This larger difference between the maximum and minimum total pressure ratios, for the numerical 
solution, produces a larger distortion value. 

The higher mass flow rate case shown in Fig. 8 b shows the same trend in the distortion value where DPCP avg = 
0.057 for the numerical simulation and 0.048 for the experiment. The difference in the low and high total pressure 
values between the CFD and experiment can also be seen by comparing these contour plots. Despite the difference in 
the distortion values, the total pressure recovery between the experiment and numerical solutions compare well. The 
low mass flow rate case in Fig. 8 a shows a Pravg/ (Pr ) °° = 0.958 from CFD results as compared to the experimental 
value of 0.957. The high mass flow rate case in Fig. 8 b shows a Pravg/ ( Pr )«, = 0.961 for the numerical simulation and 
0.963 for the experiment. 

D. Controlled: Low Mach and Low Reynolds Number Case 

The BART experiment evaluated both passive and active flow control devices. A numerical simulation of one passive 
control device configuration was compared to the experiment as well as a single active flow control configuration. 
In the active flow control configuration, several numerical simulations were compared to varying actuator mass flow 
rates. 

7. Passive Control Devices: VG Vanes 

While the focus of the BART experiments was active flow controls, a configuration using passive flow control devices 
was tested for performance comparisons and for validation of flow solver codes. The vanes were located on the same 
plane as the VG jets with eight vanes located at the bottom of the duct and six on the upper part of duct. The vanes 
had a tapered leading edge with a height of 3/16 inches and a length of 3/4 inches and were placed at an angle of 
20 deg with respect to the center line of the duct. The vanes were located along the VG plane shown in Fig. 4 and 
angled in such a way as to induce the flow away from the bottom of the duct centerline toward the top. The vanes 
used in the numerical simulation were rectangular with the same height and area. It was discovered by Allan 1 that 
when modeling trapezoidal vanes submerged in a boundary layer with rectangular vanes, that one should use the same 
height as well as the same area. A comparison of the numerical simulation of the BLI inlet with VG vanes and the 
experiment is shown in Fig. 9. This figure shows the total pressure ratio contours at the AIR The numerical solution 
clearly shows the low total pressure regions from the stream wise vortices created by the upper VG vanes. There is 
also a region on the lower side of the AIP where the stream wise vortices from the bottom VG vanes and the boundary 
layer merge. These regions can be related to the experimental contour plots as well as the interpolated CFD data. The 
numerical results predict a flow distortion of DC6 0 = 0. 134 and compare very well to the experimental value of 0. 126. 
This comparison shows how the numerical solution can be used to gain a better insight into the flow features generated 
by the flow control devices. 
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2. Active Control Devices: VG Jets 



Numerical simulations for the active flow control devices are 
compared to the best jet configuration found in the experi- 
ments of Gorton et al. The best VG jet configuration was 
found by trying different combinations of the 32 possible jet 
locations. The best configuration was defined as the config- 
uration which was able to obtain a desired DC60 value of 
0.10 using the lowest possible mass flow rate. It was found 
by Gorton et al. that the best configuration was obtained by 
using fourteen jets on the bottom of the inlet. The VG jets 
were aligned 90° to the flow and inclined 30° to the surface 
of the duct with a orifice diameter of 0.060 inches. The nu- 
merical solution for the VG jet case was modeled using 29 
grids resulting in a total of 11 million grid points. 

The experimental results presented here are for pulsed 
jets at 50 Hz with an 80% duty cycle; the numerical simu- 
lations assume steady blowing jets. From the experiments, it 
was found that there was less than 1% difference in the dis- 
tortion values between steady blowing jets and the pulsed jets 
at 50 Hz with an 80% duty cycle at the same mass flow rate. Therefore the numerical steady blowing jet results are 
compared to the experimental pulsed jet data with the assumption that the effects from the pulsed jets are equivalent 
to a steady blowing jet at the same mass flow rate. 

A mass flow sweep was performed for the VG jets at Moo = 0.15 and for a constant duct mass flow rate of 
1800 scfm. Figure 10 shows the comparison of the DC60 distortion values between the experiment and the nu- 
merical solutions. This comparison shows that the flow solver was able to capture the correct trend and was able to 
predict the point where DC60 starts to climb with increasing mass flow rate. It can also be seen that the numerical 
simulations under predicted the distortion values at lower actuator flow rates by an average of 13%. 

The surface Cp measured on the duct’s upper and lower centerline is shown in Figs. 11 and 12 for the controlled 
inlet case. Figure 1 1 compares the surface Cp for the low actuator mass flow rate case of 3 scfm where the inlet 
leading edge is located at v = 0 and the actuators at v = 6.6 inches. This plot shows that the numerical solution for 
the surface pressure matched the experiment with a slight offset. Figure 12 shows the surface Cp for a high actuator 
mass flow rate of 14 scfm. The CFD surface pressure distribution compares well with the experimental data with the 
numerical simulation capturing the low pressure region near the jets. 

The total pressure ratio contours for the actuator mass flow sweep are shown in Fig. 13. Overall, the contour plots 
for the interpolated CFD and the experiment compare very well. The location and size of the low pressure regions 


Figure 10. A comparison of the DC60 distortion values between 
the numerical and experimental results for the low free-stream 
Mach number BLI inlet case with VG jets. 




Figure 11. Upper and lower duct centerline surface Cp com- 
parison for the controlled BLI inlet with flow control jets. 
The data is from the BART experiment with a low actuator 
mass flow rate of 3 scfm, Moo = 0.15 and a duct mass flow 
rate of 1800 scfm. 


Figure 12. Upper and lower duct centerline surface Cp com- 
parison for the controlled BLI inlet with flow control jets. 
The data is from the BART experiment with a high actuator 
mass flow rate of 14 scfm, Moo = 0.15 and a duct mass flow 
rate of 1800 scfm. 
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Figure 13. A comparison of the numerical and experimental results for the BLI offset inlet, using VG jets at various mass flows, for the low 
Mach number case. These contour plots show the total pressure ratio at the AIP. The CFD results were interpolated onto the 120 probe 
locations used in the experiment. 
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Figure 14. Total pressure ratio contours of the numerical simulations for the BLI inlet with 14 control jets for varying jet mass flow ratios 
where Moo = 0.784 and Reo = 13.8 million. The location of the control jets are shown and are at the same location as the jets in the low 
Mach case where they are on the bottom of the duct at x/L = 0.34. 


on the sides of the AIP match up well with the experiment. The shape of this region does differ slightly for the low 
actuator mass flow rates. As the mass flow is increased, the shape of this low pressure region starts to differ more 
between CFD and experiment. The difference in the size and shape of this region is most likely related to the size of a 
stream wise vortex as seen in the high resolution CFD results. Finally, the total pressure contours at the bottom of the 
AIP compare very well with the experiment for each of the mass flow rates. 


E. Controlled: High Mach and High Reynolds Number Case 

The performance of the control jets were evaluated for the high Mach and Reynolds number case using the same CFD 
approach as in the low Mach number case. These simulations will provide guidance on the placement of the control 
jets for a future high Mach number BLI inlet experiment. As in the low Mach case the jets are aligned 90° to the 
flow and inclined 30° to the surface of the duct with a orifice diameter of 0.040 inches. These simulations are for the 
2.5% scale BLI inlet geometry used in the 0.3 -Meter Transonic Cryogenic Tunnel test. The performance of the jets 
are evaluated for the Moo = 0.833 case in Fig. 8(b). Of all of the flow conditions considered, this case had the largest 
DPCP avg distortion value. In all of the high Mach number controlled simulations the back pressure for the inlet is fixed 
at baseline value. Therefore, there was no attempt to adjust the exit pressure in order to maintain a constant inlet mass 
flow rate for varying jet mass flows, as was done in the low Mach and low Reynolds number cases. 

Several actuator configurations were evaluated for the high Mach number case and will be referred to as config- 
urations A through G. Jet configuration A, shown in Fig. 14, was the first to be evaluated and is the same jet layout 
used in the low Mach number cases. The actuator placement for configuration A consisted of a single row of jets at 
x/L = 0.34, where there are 7 jets on each side inlet centerline. These jets are aligned such that they are blowing from 
the centerline outwards. The simulation of the jets was performed at four different mass flow rates. The total actuator 
mass flow rates are given in terms of a jet mass flow ratio (MFR) which is defined in (4). 


Jet MFR = 


Total Jet Mass Flow Rate 
Inlet Mass Flow Rate 


( 4 ) 


The highest jet MFR for configuration A was 0.69% and demonstrated a reduction of the baseline DPCP avg distortion 
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Figure 15. Total pressure ratio contours of the numerical simulations for the BLI inlet with 14 control jets for varying jet mass flow ratios 
where Moo = 0.784 and Reu = 13.8 million. The jets are placed near the entrance of the inlet at x/L = 0.13. 


value from 0.0579 to 0.0410. The Z)C60 distortion descriptor was reduced from 0.536 to 0.367. The pressure recovery, 
Pravg / (P t) ooi was not effected by the control jets and remained at about 0.965. The low jet MFR= 0.23% shows 
a spreading out of the concentrated low pressure region seen at the bottom of the AIP in the baseline case. The 
intermediate jet MFR= 0.44% case shows the low pressure region smeared out along the bottom of the AIP with very 
little thinning of the ingested BL. The 0.52% jet MFR case starts showing a pooling of a low pressure region on the 
sides of the AIP and a slight thinning of the low pressure region. The highest MFR case of 0.69% shows larger low 
pressure regions on the side of the inlet AIP with a thinning of the low pressure region on the bottom of the inlet. 

Figure 14 also shows the ratio of the capture area, Aoo, to 
the inlet AIP area, A aip, where A^ is defined by (5). 


Aoo poo Voo = A aip p aip Vaip 


(5) 


The area averaged density and velocity at the AIP are given 
by p aip and Vaip respectively. The baseline case in Fig. 14 
has Aoo /A aip = 0.725 which shows that the stream tube ahead 
of the inlet is smaller than the inlet exit area indicating a de- 
celeration of the flow ahead of inlet. For the high jet MFR 
case, the Aoo/Aaip value of 0.719 indicates a slightly lower 
inlet mass flow rate than the baseline case. 

Actuator configuration B is shown in Fig. 15 where the 
row of 14 jets are moved upstream near the entrance of the 
inlet atv/L = 0.13. The high jet MFR case of 0.69% shows 
a much lower DC6 0 value of 0.277 as compared to config- 
uration A where DC6 0 = 0.367. Figure 16 shows the area 
averaged 60° sector total pressure ratio as the sector is ro- 
tated about the AIP for the baseline and jet configuration A 
and B cases. The ’worst’ sector average total pressure, ( Pravg)crit , for all three cases is the minimum point at 180° 
which is the bottom of the AIP. This comparison of the mean sector pressures shows how the jets in configuration B 
improve the distortion at the bottom of the inlet as compared to configuration A and the baseline. 



Figure 16. A comparison of the area averaged total pressure 
ratio for a 60° sector rotated about the AIP between the baseline 
case and the controlled cases for jet configuration A and B. 
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Figure 17. Total pressure ratio contours of the numerical simulations for the BLI inlet with 28 control jets for varying jet mass flow 
ratios where Moo = 0.784 and Reu = 13.8 million. This case has two rows of jets on the bottom surface of the BLI inlet at x/L = 0.13 and 
x/L = 0.34. 


A comparison of total pressure ratio contour plots for configuration A and B shows how the jet actuators placed 
closer to the entrance of the inlet improved the average pressure at the bottom of the inlet. However, configuration B 
is less effective in moving the low pressure along the sides of the inlet. It also shows how the DPCP avg only varied 
slightly (less than 5%) where DC60 was reduced by about 30% when moving the jets closer to the inlet entrance. 

Configuration C is a combination of the two jet actuator locations, and has 28 jet actuators as shown in Fig. 17. The 
low jet MFR case is 0.63% with a DC60 = 0.429 and a DCPC avg = 0.0553. These distortion values are much higher 
than the two previous cases at the same jet MFR and with half the number of jets. This indicates that the momentum 
of the jets and not the total mass flow rate is the most important parameter for the effectiveness of the control jets. 
Increasing the jet MFR to 1.38% for this case results in DC6 0 = 0.245 which is 15% lower than the single jet row 
case in Fig. 15 but at twice the jet MFR. The DPCP avg for the two jet row case at the high jet MFR in configuration 
C is 29% lower than high jet MFR for configuration B. This shows how improving one distortion descriptor does not 
necessarily result in an improvement of the other descriptor. 

Figure 18 shows the location of the jets and numerical results for configuration D where 28 jets are located on the 
bottom and side of the inlet at x/L = 0.13. The 1.37% jet MFR case had a DC60 = 0.255 and a DPCP avg = 0.0341. 
Compared to jet configuration C, this configuration has higher distortion values for the same amount of jet MFR. This 
shows that adding the jets on the bottom of the inlet downstream is more effective in reducing distortion than adding 
them on the side of the inlet upstream. Since the ’worst’ 60° sector is on the bottom of the inlet, it makes sense that 
adding jets on the side would be less effective than putting them on the bottom. 

Configuration E is identical to D but with a downstream row of jets added on the duct bottom for a total of 42 
jets, as shown in Fig. 19. The low jet MFR case of 1.34% had a DC6 0 = 0.334 and a DPCP avg = 0.0391 which are 
much higher than configuration C for the same jet MFR. The high jet MFR for configuration E was 2.08% producing 
DC6 0 = 0.224 which is only 9% lower than configuration C at a jet MFR of 1.38%. Configuration E for the MFR 
case of 1.34% had a DPCP avg = 0.0391 which is also only 9% lower than configuration C at 1.38% jet MFR. This 
comparison between configurations C and D shows how adding the side jets only results in a small improvement of 
the flow distortion while significantly increasing in the total jet MFR. 

Jet configuration F shown in Fig. 20 has a more distributed pattern of jets on the bottom and side of the inlet. This 
case has 30 jets on the bottom and 10 jets on the side of the inlet near the comers. The total pressure contour plots 
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Figure 18. Total pressure ratio contours of the numerical simulations for the BLI inlet with 28 control jets for varying jet mass flow ratios 
where Moo = 0.784 and Reu = 13.8 million. This case has one row of jets near the entrance at x/L = 0.13 with 14 jets on the bottom and 14 
jets on the side of the BLI inlet. 
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Figure 19. Total pressure ratio contours of the numerical simulations for the BLI inlet with 42 control jets for varying jet mass flow 
ratios where Moo = 0.784 and Reo = 13.8 million. This case has two rows of jets on the bottom surface of the BLI inlet at x/L = 0.13 and 
x/L = 0.34 with a row of jets on the sides near the entrance at x/L = 0.13. 
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Figure 20. Total pressure ratio contours of the numerical simulations for the BLI inlet with 40 control jets for varying jet mass flow ratios 
where Moo = 0.784 and Reo = 13.8 million. This case has five rows of jets with each row having six jets on the bottom and two on the side 
of the BLI inlet. 
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Figure 21. Total pressure ratio contours of the numerical simulations for the BLI inlet with 56 control jets for varying jet mass flow ratios 
where Moo = 0.784 and Reo = 13.8 million. This case has eight rows of jets on the bottom surface of the BLI inlet with the first row having 
14 jets and the other rows with six jets concentrated at the center of the BLI inlet. 
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Figure 22. Summary of DC60 versus the total jet MFR for all of 
the jet configurations for the high Mach and Reynolds number 
cases. 



Figure 23. Summary of DPCP avg versus the total jet MFR for all 
of the jet configurations for the high Mach and Reynolds num- 
ber cases. 


in Fig. 20 show that this actuator pattern produces two large low pressure regions on the side of the inlet AIR The 
distortion for the high jet MFR case shows that DC60 is only slightly higher for the same jet MFR in configuration E. 
Likewise, DPCP avg is 14% higher when compared to the same jet MFR in configuration E. Jet configuration F turns 
out not to perform as well as configuration E because of the higher distortion values and the larger low pressure regions 
on the side of the inlet. 

Configuration G looked at concentrating the jets at the bottom center of the inlet where the distortion is the highest. 
Figure 21 shows the location of the eight rows of jets for this configuration. This configuration had a total of 56 jets 
where the first row had 14 jets and the next seven rows had six jets concentrated near the bottom centerline of the 
inlet. The low jet MFR case was 1.82% and had a DC60 distortion of 0.226 and a DPCP avg of 0.0330. The jets on 
the centerline do a good job of clearing out the center but create large low pressure regions on the side of the inlet. 
These regions produce the worst 60° sector. Increasing the jet MFR to 2.32% improves the distortion at the bottom 
of the AIP and moves the low pressure regions higher along the sides. The 60° average sector pressure for these 
low pressure regions improves slightly lowering DC6 0 to 0.211. The 2.32% jet MFR case also reduces the DPCP avg 
to 0.0302. Increasing the jet MFR to 2.87% increases the size of the low pressure regions on the side of the inlet 
without changing the flow on the bottom of the AIP. This increase in the size of the low pressure regions increases the 
distortion values where DC6 0 = 0.246% and DPCP avg = 0.0343. All three controlled cases caused a small reduction 
in the pressure recovery compared to the baseline case. 

A plot of the distortion descriptors versus the total jet MFR for each of the jet configurations is shown in Fig. 22 
and 23. The comparison of the DC60 distortion in Fig. 22 shows how none of the configurations were able to reduce 
DC60 below 0.2. Comparing configuration A with B indicates that moving the row of jets from x/L = 0.34 to 0.13 
greatly improved DC6 0 for the same number of jets at the same jet MFR. The other configurations indicate, that by 
adding more jets and increasing the total jet MFR, predicted DC60 can be decreased another 0.065 but at the cost of 
tripling the total jet MFR. The lowest predicted DC6 0 was 0.21 1 for jet configuration G with a total jet MFR of 2.32%. 
The best jet configuration for a jet MFR under 1% was configuration B where DC 6 0 = 0.277 for a jet MFR= 0.69%. 

Figure 23 shows the trend of DPCP avg versus jet MFR for the controlled high Mach number cases. This plot shows 
that DPCP avg bottoms out at 0.0281 for a jet MFR of 2% using configuration E. This figure also indicates that there is 
little difference between DPCP avg for configurations A and B and how the combination of the two, in configuration C, 
continues the decreasing trend of DPCP avg for increasing jet MFR. Figure 23 also shows how the low jet MFR cases 
for configuration C and D did not perform as well as A and B for a given jet MFR. 

A summary of the distortion descriptors in terms of total jet momentum coefficient, Cm jet, are shown in Figs. 24 
and 25, where the total jet momentum coefficient is defined in (6). 
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Figure 24. Summary of DC60 versus the total jet momentum 
coefficient, C^jet? for all of the jet configurations for the high 
Mach and Reynolds number cases. 



Figure 25. Summary of DPCP avg versus the total jet momentum 
coefficient, C^jet? for all of the jet configurations for the high 
Mach and Reynolds number cases. 
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These figures indicate that distortion is more a function of the total jet momentum as opposed to the jet mass flow 
ratio. As shown in Fig. 24, all but two configurations tend to fall on the same curve. Configurations B and G have 
lower distortion values for the same jet momentum as compared to the other configurations. Figure 25 plots DPCP avg 
as a function of the total jet momentum and all of the cases tend to fall on the same curve. Only the low jet MFR case 
for configuration D seems to fall off of this curve with higher DPCP avg values. These figures indicate a trend where 
the distortion starts to flatten out for increasing total jet momentum. 


VI. Summary 

The objective of this research was to evaluate the ability of a RANS flow solver to predict the flow quality at the 
location of the engine face for a BLI offset inlet for baseline and flow controlled cases. Once the numerical approach 
was validated, the performance of VG jets, for the BLI inlet at high Mach and Reynolds numbers were evaluated. These 
numerical simulations will then provide guidance on the design of a future wind tunnel test. The second objective of 
the study was to provide insight into the flow physics of the flow control actuators for the development of a source 
term actuator model. 

The BLI offset inlet geometry was generic in nature so that it would be an open geometry and still be representative 
of an inlet for a commercial version of the BWB aircraft. Numerical simulations were compared to two different 
experimental tests using this inlet geometry. The experiments were setup such that the inlet was ingesting a boundary 
layer with a height of at least 30% of the inlet height. The experimental data consisted of a high Mach and Reynolds 
number test for a baseline inlet (i.e. no flow control devices) performed in NASA Langley’s 0.3-Meter Transonic 
Cryogenic Tunnel. A low speed test in NASA Langley’s Basic Aerodynamics Research Tunnel (BART) provided data 
for the BLI inlet with passive and active flow control devices. 

Overall, the numerical results compared well with the experiments for both the baseline and controlled cases. 
The baseline flow distortion for the high Mach and Reynolds cases were over predicted by the flow solver. This was 
related to the flow solver predicting a slightly lower minimum total pressure than was measured in the experiment. 
Despite this difference in the distortion, the flow solver was able to predict the pressure recovery well. Overall, the 
flow features in the contour plots of the total pressure compared well with experiments. 

The low Mach and Reynolds number comparisons showed that the flow solver under predicted the baseline flow 
distortion by 13%. However, the contour plots showed very similar flow features and compared very well with the 
experiments. A comparison of the duct flow using passive flow control devices, VG vanes, showed good agreement 
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between the CFD and the low Mach number experimental data. Simulations of the VG jets were compared to an 
actuator mass flow sweep performed in the low Mach number experiment. This comparison indicated that the flow 
solver was able to predict the inlet distortion for varying actuator mass flow. The flow solver was also able to predict 
the actuator mass flow rate at the minimum distortion point. 

The performance of the VG jets for the BLI inlet at high Mach and high Reynolds numbers were evaluated for 
several different actuator locations. These simulations indicated that a minimum DC60 of 0.211 could be achieved at 
a cost of 2.32% of the inlet mass flow. The simulations also predicted a minimum DPCP avg of 0.0281 using 2.08% of 
the inlet mass flow. From an inlet design point of view, a DC60 under 0. 10 is desirable while using less than 1% of the 
inlet mass flow. The simulations for the high Mach number case indicate that more research needs to be performed on 
the location of the actuators in order to reduce the flow distortion further. Also the simulations showed that an inlet 
mass flow of 2% or more may be needed to acieve a DC6 0 value less than 0.10 

The cost of fully griding the jet actuators was very high so only a few actuator patterns could be explored in this 
study. The results from these simulations will be used to validate a source term modeling approach which will greatly 
reduce the cost of modeling the jet actuators for the BLI inlet. 
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